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FLUCTUATIONS IN THE HOMOGENIZATION OF SEMILINEAR 
EQUATIONS WITH RANDOM POTENTIALS 
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Abstract. We study the stochastic homogenization and obtain a random fluctuation theory for 
semilinear elliptic equations with a rapidly varying random potential. To first order, the effective 
potential is the average potential and the nonlinearity is not affected by the randomness. We then 
study the limiting distribution of the properly scaled homogenization error (random fluctuations) 
in the space of square integrable functions, and prove that the limit is a Gaussian distribution 
characterized by the homogenized solution, the Green’s function of the linearized equation around 
the homogenized solution, and by the integral of the correlation function of the random potential. 
These results enlarge the scope of the framework that we have developed for linear equations to the 
class of semilinear equations. 
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1. Introduction 

1.1. Motivation and overview. We study the asymptotic stochasticity of solutions to the semi- 
linear elliptic equation with random potential 

{ — Au^ + Qsix, uj)u‘^ + = 9{x), X G D, 

u^ = 0, xedD, 


and characterize the limiting distribution of the random fluctuations (correction to homogeniza¬ 
tion). Here, D is an open bounded domain in d = 2,3, with smooth (say boundary dD. 
The (linear) potential qe{x,uj) is highly oscillatory and is of the form q{j,uj), where 0 < e <C 1 
denotes the scale of oscillations. We assume that q{x,uj), the potential function before scaling, is 
a stationary ergodic random field on a probability space (H, ju^p). xhe nonlinearity / is assumed 
to satisfy conditions that, essentially, guarantee that for each a; G H and g G the above 

equation admits a unique weak solution. 

Under mild conditions, such as, e.g., stationarity and ergodicity, on the random potential q{x,u) 
and under natural structural conditions on the nonlinearity f{u), we prove in Theorem 13.31 below 
that the above semilinear problem homogenizes as e —?■ 0 and that the effective potential is in fact 
given by averaging and the nonlinearity term is not changed. In other words, converges strongly 
in L^{D) and weakly in H^{D), for a.e. a; G H, to the solution u of the effective equation 

{ — An + qu + f(u) = gix), x G D, 

( 1 . 2 ) 

n = 0, X G dD, 

where q = ¥.q is the average of q{x,uj). The homogenized equation hence is deterministic and has 
smooth (non-oscillatory) coefficients. 
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The main contributions of this paper is an estimate of the size of the homogenization error — u, 
say in the Lp‘[D)) norm, and more importantly, a study of the law (probability distribution) 

of the random fluctuations — n, viewing them as a random element in the Hilbert space L?‘{D). 
Such asymptotic estimates of stochasticity beyond the homogenization limit find applications in 
uncertainty quantification and Bayesian formulation of PDE-based inverse problems; see e.g. [HIES]. 

The study of the limiting distribution of the homogenization error goes back to m, where 
the Laplace operator with a random potential formed by Poisson bumps was considered. General 
random potential with short range correlations was considered recently in [T], and in [U [5] for 
other non-oscillatory differential operators with random potential. Long range correlated random 
potential was considered in [2]. When random elliptic differential operators are considered, the 
limiting distribution of the homogenization error was obtained in [9j for short-range correlated 
elliptic coefficients, and in [3] in the long-range correlated case; all in the one dimensional setting; 
we refer to [la El EH EH EO] for important recent advances in this direction in higher spatial 
dimensions under strong structural assumptions on the probability distribution of the coefficients. 

The main results of this paper show that the general framework developed by the authors in 
[HElElllB], for linear equations with random potential, essentially applies to the semilinear equation 
dni). Let us briefly explain the main ideas. In the linear framework, e.g., for (jl.ip with / = 0, 
let us denote the fundamental solution operator (Green’s operator) of the Dirichlet problem for 
— A -|-g by Q. Then the homogenization error — u admits the following expansion formula: 

— U = —Qv^U + “ u). (1.3) 

Here, — q denotes the random fluctuation of the potential. We also have — u = —G^VeU, 

where G^ is the solution operator associated to —A + q^. As long as q^ > —Ai, the first eigenvalue of 
the Laplace operator, we obtain a bound on — u in L^(H, Lp‘{D)) provided that we can estimate 
the correlation function of The leading term in the right-hand side of (|1.3p is the first one. 
The last term is smaller because both — u and the operator Gi^eG are small. The middle term is 
also of smaller order but for a different reason: one can verify that Gi^eGi'eU ~ '^{G'l^eGveu) is small, 
provided we can estimate fourth order moments of and ¥.{Gi^eGveu) is also small in dimensions 
d = 2,3. 

For the nonlinear equation, ()1.3p has to be modified. In fact, as was already pointed out in [12] 
without technical details, the leading term oi — u in distribution should be given by —Gu^eU-, 
where Gu is the solution operator associated to the linearization of the homogenized equation around 
the homogenized solution u. In this paper, we derive an expansion formula for — u (see dO]) 
below), which plays the role of (II.3p . In fact, G above is replaced by Gu and there are two more 
terms involving the nonlinearity /. To overcome the difficulty caused by the nonlinearity, we first 
establish uniform and L°° estimates on solutions to (HI) that are independent of e and uj. These 
uniform bounds, together with regularity assumptions on /, allow us to estimate f{u^) — f{u) as 
— u and estimate f{u^) — f{u) — f'{u){u^ — u) as (u^ —tt)^. As a result, we show that the nonlinear 
terms in the expansion oi — u do not contribute to the limiting distribution. 

Finally, the picture of the size and the limiting distribution of the homogenization error — u 
is as follows; see Theorem 12.21 and Theorem 16.21 below. In dimension d = 2,3, for short range 
correlated random potential, — n is of order ez in the L^(H,L^(L>)) norm, and the distribution 
in the space L^{D) of the normalized error e“ 2 (u^ —tt) converges weakly to a Gaussian distribution 
on LP‘{D), with correlation kernel determined by the homogenized solution u, the Green’s function 
associated to Gu, and the integral of the correlation function of v{x,uj). In the long range case, for 
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potentials that are functions of long range correlated Gaussian random fields, the size of — u 
is of order where a E (0,d) is the decaying rate of the correlation of the underlying Gaussian 
field. The limiting distribution of [vF — rt) is a long range correlated Gaussian, with correlation 
kernel depending also on a and certain parameters in the model of the random potential. 

The rest of this paper is organized as follows: in section [2] we make precise the main assumptions 
on the random potential and the nonlinearity in (II.ip . and state the main theorems. In section [3l 
we establish well-posedness of and more importantly, uniform and L°° bounds that are 
independent of e and co. These results facilitate the homogenization proof. Section 0] and [5] are 
devoted to the proof of the main results, where we characterize how the homogenization error scales 
in energy norm, and determine the limiting distribution of the scaled error. We recall the controls 
on the terms in the expansion formula that are the same as in the linear setting, without detailing 
the proof, and concentrate on how to apply the uniform L°° bound to deal with the nonlinear 
terms, which are new. Finally, we make some comments in section [6] and discuss some directions 
of further studies. 


2. Assumptions and Main Results 

2.1. Assumptions. We first state the main assumptions on the nonlinearity / and the random 
potential q{x,uj). Let Xi{—A;D), which is henceforth simplified to Ai(—A) or even Ai, be the first 
eigenvalue of the Dirichlet Laplacian operator on D, that is, 

Ai(—A) := inf |y \Dw\^ dx : w G Hq{D), J = l| . ( 2 - 1 ) 

Since D is bounded, we note that Ai(—A) > 0. The first set of main assumptions on q and / are: 

[Al] q{x,u}), X E M'^, is a stationary ergodic random field on (R, P). This means, there exist 

a family of P-preserving ergodic group of transformations {t^ : 12 —)• i.e. 

¥{txA) = P(A) for all x E and A G F. 

TyA = A for all y E implies that P(A) E {0,1}, 

and a random variable q such that q{x,Lo) = q{TxOj). 

[A2] There exist M >1 such that 

|g(a;)| < M, for all a; E 12. (2-2) 

Moreover, there exists 7 > 1 such that |/(s)| < C{\ + Isp) for all s E M, and if d = 3, we 
further assume that 7 < 2d/(d — 2) — 1 = 5. 

[A3] / : M —>■ M is continuously twice differentiable, and for some c > 0 , 

Al + q{u}) + f'{s) > c, for all a; E 12 and s E M. (2.3) 

The assumptions [A2] and [A3] guarantee that for each cj E 12, the heterogeneous semilinear 
equation o admits a unique weak solution. The stationarity and ergodicity of q{x,uj) further 
makes sure that this problem homogenizes, in the limit e —>• 0 , to the effective equation (| 1 . 2 p . 
Notice that [A2][A3] still hold if we replace ^by q; hence, the homogenized problem is also uniquely 
solvable. 

To get quantitative estimates on the homogenization error F — u, and to find the limiting 
distribution of this error after proper scaling, we need more information on the random potential 
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q{x,io), or equivalently the fluctuation 

uj) := q{x, ui) —Kq = q{x, uj) — q. (2.4) 

In most of the paper, we assume that q{x,uj) is a short-range correlated random field satisfying 
certain fourth-order moment estimates. Let R{x) = ¥.{y{x,uj)v{d,uj)) = E(i^(x + be 

the auto-correlation function of the stationary field u. The term “short range correlation” refers to 
the condition that the correlation function R{x) is integrable, i.e. 

cr^ := / R{x)dx < oo. (2-5) 

jRd 

Note that by definition, R{x) is a nonnegative definite function in the sense that for any positive 
integer n, for any n-tuples (xi, • • • , Xn), the matrix formed by {R{xi — Xj))ij=i^... is nonnegative 
deflnite. By Bochner’s theorem, the integral of R, i.e. cr^, is nonnegative. Throughout the paper, 
we assume also that > 0. 

We will need an estimate for (mixed) fourth order moments of u later. To simplify the presenta¬ 
tion, we impose a stronger condition using the notion of “maximal correlation coefficient”, which 
quantifies how fast the correlation of u decays. Let C the set of compact sets in and for two 
sets Ki, K 2 in C, the distance d{Ki, K 2 ) is defined to be 

d{Ki,K 2 )= min \x - y\. 
x£Ki,y£K 2 

Given any compact set iL C C, we denote by Rk the cr-algebra generated by the random variables 
{i^{x) : X G K}. The maximal correlation coefficient ^ of is defined as follows: for each r > 0, 
Q{r) is the smallest value such that the bound 

E((/?i(z/)(/?2(i^)) < £'(r)yE7^f(^^)yE7^2^ (2.6) 

holds for any two compact sets Ki,K 2 G C such that d{Ki,K 2 ) > r and for any two random 
variables of the form i = 1,2, such that is J-Xi-measurable and Kipi{u) = 0. 

The further assumption we have on q[x,uj) is: 

[S] The maximal correlation function satisfies G L^(]R+, r'^“^dr), that is 

g2 {r)r‘^~^dr < 00 . 

Assumptions on the mixing coefficient g of random media have been used in [I10II6]; we refer 
to these papers for explicit examples of random fields satisfying the assumptions. Note that the 
correlation function R{x) of v can be bounded by g. For any x G 

|i?(x)| = \Ei'{x,uj)iy{0,uj)\ < ^(|x|)Var(g). 

By [A2][A3], q, and hence its variance are bounded. Since we may assume g G [0,1] (hence g < s/g), 
[S] implies R G ^^(M'^) and q{x,uj) has short range correlations. In fact, |i?| 2 G so roughly 

speaking, |i?(x)| < I/|x|“ with a > 2d. Therefore, [S] is a stronger condition than q being short 
range correlated, which would just require the decay rate a > d. This faster decay of correlation 
is not necessary for the main results of this paper to hold (see Remark 12.41 below). But using the 
assumption [S], we can simplify significantly the following fourth order moment estimates of the 
random potential u{x,uj). 
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For any four points x, y, t and s in define 


^uix,y,t,s) :=Eu{x)i^{y)iy{t)u{s) - (Eiy{x)u{y))(Eiy(t)i^{s)). 


(2.7) 


Were ly a Gaussian random field, its fourth order moments would decompose as a sum of products 
of pairs of R and the above quantity simplifies to a sum of two products of correlation functions. 
This property does not hold for general random fields. However, we have the similar estimate: 


\^„{x,t,y,s)\ < t?(|x - tD'dily - s|) - s|)i9(|y - t\) 


( 2 . 8 ) 


where 7?(r) = {KQ{r/3))^, K = 4:\\iy\\L°°{nxD)- We refer to [TB] for the proof of this lemma. 
Estimates of this type based on mixing property already appeared in [T]. 


Notations. We simplify the notation of Lebesgue spaces U’{D) on D to when this is not 
confusing; in particular, if Lebesgue space on the probability space H is concerned, we make the 
dependence explicit. We use for the Sobolev space and we use 77®, s G (0,1), for the 

fractional Sobolev space which is the closure of C^{D) in the norm 


\u\ 


H^(K) 


:= m 


L^{K) 


+ 


/ 


\u{x) - n(y)d 


\x ■ 


_ ^|c/+2s 


■dxdy. 


See m for more reference on 77®. Throughout the paper, C denotes various bounding constants 
that may change line after line and we say C is universal when it depends only on the parameters 
in the main assumptions above. For the functions q, v and 7?, we use subscript e to denote the 
scaled function, as in q^ix) = q (|). 

2.2. Main results. The first main theorem concerns how the homogenization error scales. 

Theorem 2.1. Let D C be an open bounded domain with boundary dD, n® and u be the 
solutions to dni and m respectively. Suppose that [Al][A2][A3] and [S] hold, g G L‘^{D) and 
d = 2,3. Then, there exists positive constant C, depending only on the universal parameters and 
||( 5 :||j;^ 2 , such that 

¥.\\u^-u\\\2 <C{\\g\\L2)e‘^. (2.9) 

This theorem provides LP‘{TI,LP‘{D)) estimate of tt® — u, and its proof is detailed in section [H 
So roughly speaking, the size of the homogenization error is of order ez. We hence consider the 
limiting law of the rescaled random fluctuation — u), and prove: 

Theorem 2.2. Suppose that the assumptions in Theorem \2 .71 hold. Let a be defined as in (USD 
and Gu{x,y) be the Green’s function defined by the linearized equation around u; see (|3.6p below. 
Let W(y) denote the standard d-parameter Wiener process. Then as e ^ 0, 


U U distribution 


(7 


[ Gu{x,y)u{y)dW{y), in L^{D). 
J D 


( 2 . 10 ) 


The proof can be found towards the end of section [5l 


Remark 2.3. The integral on the right hand side of (I2.10p is understood, for each fixed x, as 
Wiener integral in y with respect to the multiparameter Wiener process W{y). Let X denote this 
integral. For d = 2,3, because the Green’s function Gu{x,y) is square integrable (see below), X 
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is a random element in L'^(D). In particular, for any (p E L‘^{D), the inner product {p>,X) has 
precisely the Gaussian distribution AA(0, cr^) with mean zero and variance 

■= f {Qu^fv?dy. ( 2 . 11 ) 

Jd 

Remark 2.4. We make some comments on the main assumptions of this paper. In view of the 
variational formulation of assumption [A3] is natural because it guarantees that the mini¬ 

mizing functional is coercive and convex, as we show in the next section. The full zero order term 
f^{u) := QeU -|- f{u) can be viewed as a reaction-diffusion nonlinearity. As an example, when / is 
of bistable type, say f{u) = u{u — l)(u — 9) for some 9 E (0,1), and q is large enough, all of the 
requirements in [A3] are satisfied. 

As pointed out earlier, the assumption [S] imposes stronger decay rate on the correlation function 
R than it being integrable. This stronger condition is mainly imposed to simplify the control of the 
fourth order moment in (|2.8p . We refer to [5] for an alternative way to control terms like 
which allows decay rate |R(x)| < |x|“", a > d, at the cost that there are more pairs of products of 
controlling functions on the right hand side of ( 12 . 8 p . 

We state and prove the main theorems assuming that the random potential q{x,uj) has short 
range correlation. Nevertheless, our approach also works for some long range correlation setting; 
see the last section of this paper. 

3. Homogenization and Uniform Estimates 

3.1. Well-posedness and uniform estimates. We record here first the well-posedness result for 
the heterogeneous semilinear equation O)- 

Lemma 3 . 1 . Under the assumptions [A2][A3], for each fixed lo £ Q, and e E (0, 1), g E , there 
exists a unique weak solution in Hq{D) to (jl.ip . Moreover, there exists C > 0 independent of e 
and oj, and 

\\u"\\h^ <C{\\g\\H-i)- (3.1) 

This result is classical; nonetheless, we briefly outline the proof for the sake of completeness. 

Proof. For fixed ui and s, one can recast (jl.ip as the Euler-Lagrange equation associated to the 
variation problem of minimizing the following (nonlinear) energy functional 

■= + -^9 dx, (3.2) 

over the set Hq{D), where F is an antiderivative of /, e.g. F{v) = Jq f{s)ds. 

The assumption |/(m)| < C'(l-l-lup) implies that |E(u)| < C'(l-|-|u|'>'’''^). By Sobolev embedding, 
if d = 2, H^{D) is embedded in LP{D) for any p E [l,oo) and if d = 3, H^{D) is embedded in 
L^{D) which is included in since 7 < 5. In both cases, F{v) is integrable. So, the integral 

above is well defined on Hq{D). 

The assumption Ai -|- g-|- f'{s) > c > 0 for all s E M further implies that U, as a functional on 
Hq{D), is strictly convex. Moreover, we can choose a > 0 small so that {I — a)Xi + q+ f {s) > § > 0 
still holds for all s and cu. Integrate in s on both sides of this inequality; we get 

i ((1 - a)Ai qe{x,u)) -h F{s) > /(0)s. 
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Using this fact and the Poincare inequality ||Vu ||^2 > -^ 11 ^ 11^25 we get 


1 


1 


a , 


ID 


-\Vv\ +-qe{x,uj)v + F{v)dx> / —\Vv\ +-\v\ dx+ / f{0)vdx 


> 


L 


D 2 
a 


Id 


|Vup + -|up dx — 


2 |/( 0 )|^ 


(3.3) 


\D\ 


/£) 2 8 c 

where \D\ is the volume of the domain D. This shows that F is also bounded from below. As a 
result, (|3.2p and hence admits a unique solution for each e and lo. 

For the uniform bound, we solve, any g E the deterministic Poisson problem 


-Aug = 2g, in D, 

with boundary condition Ug = 0 on dD. Clearly we have ||%||hi < for some universal 

constant C. The weak formulation of Ug then yields, 

= J (^ 3 )^ +^(ug) - ffUg dx = ^ J^qe{x,uj)uldx + J^F{ug)dx. 

Now we compare F[ug] and F[u^]. Note that p.3p implies 

F[u^]>-a f [ gu^ — C. 

2 Jd Jd 

Using F[u^] < F[ug] and then Holder’s and Young’s inequalities, we have 

+ j^d^(^9)dx+ 4:C{c',Xi)\\g\\jj-i+ ^a + C, 

for some big constant C. By [A2], the integral of q^u^g is bounded by C||ttg||^ 2 ) which is further 
dominated by C'||( 7 |||^_i. In view of the growth condition on F and the Sobolev embedding as 
before, we get \\F{ug)\\Li < C{\\ug\\Hi) = C{\\g\\H-i) for some C depending also on the parameters 
in the main assumptions but is uniform in e and oj. The desired result then follows. □ 


Next, we show that the solution E L°°{D), and the bound on ||u'^||£,oo can be made independent 
of e or uj. 


Lemma 3.2. Assume that [A2][A3] hold. For any fixed s E (0,1) and u G Q, let be the 

unique solution of (ini). Then for any g E Lfi{D), there exists C > 0 independent of e and oj such 
that 

Wu^Wl^ <C{c,M,\\g\\L 2 ). (3.4) 

Proof. We use the following results from elliptic regularity theory: 

(i) Suppose that v solves the linear equation 

—Av = h in D, 

with Dirichlet boundary condition u = 0 on dD, h G with p > max{l, and v G F for some 
r E [ 1 , 00 ), then it holds that ||u||loo < Cdl/illiP + ||u||Lr-) with C also depending on r, p and D. 
We refer to m Theorem 4.2.2] for the proof of this result. 

(ii) Suppose v satisfies —Av + v = h in D with Dirichlet boundary condition on dD, h G H~^{D) 
and h G LP(D) for some p G (1, |) and d = 3. Then there exists C(p,d,D) such that ||u|| dp < 

C\\h\\Lr. We refer to m Theorem 4.2.3] for the proof of this result. 

To prove Lemma 13.21 first consider the simple situation of d = 2 and 7 > 1 arbitrary in [A2], or 
d = 3 and 1 < 7 < 4/(d — 2) = 4. For notational simplicity, we write u for . Recast the semilinear 
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equation (HH) into the form of —Au = h with h = g — f{u) — qe{x,u)u. Now assumption [A 2 ] 
imposes that ||g£tf||LP < CUttllLP uniformly in p, e and iv. If d = 2 , then by Sobolev embedding and 
the growth condition on /, we verify that h & for all p G [l,oo). If d = 3 and 7 < 4, then by 
Sobolev embedding, u ^ U with r = and r > ^. Then f{u) G L~i and ^ f • Note also that 
g G L^, QeU G and 2 > | for d = 2,3; we hence verify that /i G for p > Applying regularity 
theory (i), we conclude that 


||w||l°° < C{c,M, II^IIlz, ||w||hi) < C{c,M, \\g\\L 2 ). 

For the more general case, d = 3 and ^ (i.e. 7 G (4,5)), we need a bootstrap argument. 

Here, we mimic the proof of Theorem 4.4.1 in Cazenave m- The semilinear equation (jl.ip can 
also be recast as —Au + u = b with b{x) = g{x) + (1 — qe{x,uj))u — f{u). Again the function g and 
(1 — qs)u are in with p > | with uniform bounds, so we only need to take care of the nonlinear 
term f{u). To start, we set r = and we have u G and we note 


dy 

1 <r <—, 


2r 4 d + 2 

— =7 -< - 

d ' d -2 d -2 


4 

d -2 


= 1 , 


where we used 7 < as in [A2]. In particular, 9 := > 1. Since u G U and \f{u)\ < 

(7(1 + |tt|'^), we check that b G L~> and ^ G (1, 2 )- By elliptic regularity (ii), we obtain 

dr 

" d 7 - 2 r d - ^ 

7 

If dr > ^, we are back to the simple situation. If not, we repeat the argument using elliptic 
regularity (ii) k more times until 6^r < ^ < O^^^r. We then get u G with d^+^r > ^ and 

hence return to the simple situation, and we can conclude. We further verify that the bounds are 
uniform in e and ui in view of ()3.ip . □ 


3.2. Homogenization theory. Because the random coefficients appear only in the zeroth order 
linear term, i.e. the potential term, the homogenization theory for the equations is relatively 
straightforward. For the sake of completeness, we present the details here. 


Theorem 3.3. Under assumptions [Al][A2] and [A3], there exists an event G iF with full 
probability measure, sueh that for all uj G Hi, converges strongly in L‘^{D) and weakly in 

Hq{D), to the deterministic function u G Hq{D) that solves (|1.2n 


Proof. Owing to Lemma l3.II and Lemma 13.21 we note there exists C > 0 uniform in e and a;, such 
that \\u^\\h^{d) + ^ C. Also, by ergodic theorem (see. e.g. [TTl section 7.1]), there 

exists Hi G with P(Hi) = 1 such that q^{x,Ljj) converges weakly in to q for all ui G Hi. 

As a consequence, for each fixed a; G Hi, we can find a function u{-,io) G Hq{D) and extract a 
subsequence e(cu) — 0 , along which we have 






Q 



\e / 

where ^ and —>• denotes, respectively, weak and strong convergence, 
bound on F and the regularity assumption / G (7^(M), 


Owing to the uniform in e 


\f{u^) - f{u)\ < - u 








Here, Ci is the Lipchitz constant of / inside the domain [—2(7, 2C] and C is the uniform bound in 
(j3.4p . This implies that f{u^) f{u) in Lp‘[D). As a result, pass to the limit e —)• 0 in the weak 
formulation 




ID 


■ Vip + q — gip dx = 0, for all p E (7“(D), 


we get 


ID 


Vu ■ Vip + qup + f{u)ip — gp dx = 0, for all p E C^{D). 


This shows that the function u{-,uj) solves the equation (11.211 . Note that this problem has a unique 
deterministic solution, and hence u{x,u}) = u{x) is independent of u. As a result, for all w E Hi 
and along the full sequence e ^ 0, the convergence ^ u hold. This completes the proof of the 
theorem. □ 


3.3. Green’s function estimates for the linearized equation. Let u be the homogenized 
solution. The linearized differential operator around u, of the homogenized semilinear operator 
Cu = —Am + qu + f{u) is given by 

Cu{v) =-Av + {q + f'{u)) V (3.5) 

The Green’s function Gu{x,y), x,y G D and x ^ y, associated to £„ satisfies 

( - AGu{x]y) + {q +f'{u{x))) Gu{x-,y) = 6y, for x E L>, 

Gu{x; y) = 0, for x E dD. 

We have the following estimates: 


Lemma 3.4. Assume [A2][A3], and assume that D C is an open bounded domain with G^ 
boundary. Then there exists G > 0 such that, 


\Gu{x,y)\ < 


G 

|x — y\^~'^ 
(7(|log|x-?/|| + 1) 


for 4 = 3 
for 4 = 2 


and 


VxGu{x,y)\ < 


G 

X — y\‘^~^ 


(3.7) 


We note that by (13.41) . the potential function q + f'{u) E L°^{D). Moreover, [A3] guarantees that 
the problem remains elliptic. The first bound in (|3.7I) is immediate and the second one follows, say, 
from standard Holder regularity for gradients. Note that this is the place where regularity of dD 
is used. 


4. Quantitative Estimates of the Homogenization Error 

In this section, we determine the convergence rate of m^ —)• m in L^(H, L^(D)). Dehne := u^ — u. 
Then, it satisfies 

-Aff + qsff + fiu^) - f{u) = -i'eix,uj)u. 

Formally, the nonlinear term f{u^) — f{u) is approximated by f'{u){u^ — u). So we expect that the 
leading term in ff is given by the solution of the following equation: 

f C.uX^ = + qx^ + f'{u)x^ = in D, 

I X^ = 0) on dD. 


(4.1) 




Let Qu to be the inverse operator, i.e. the fundamental solution operator of the Dirichlet problem, 
and Guix,y) is the Green’s function. In view of [A3], these notations are well defined, and we have 

= -QuVeU. (4.2) 

Our goal is to estimate the quantity First, an estimate for L‘^{D)) is 

easily obtained from the linear theory. Next, by applying the mean value theorem to the nonlinear 
term in the equation of we show that the remainder := — x^ satisfies a linear equation 

with coefficients that depend on and u, and with on the right hand side. Therefore, by linear 
theory again, we finally obtain estimates for and hence for To use the linear theory for z^, 
however, we need the uniform (in e and ui) bound on because the coefficient of the equation for 
z^ depends on u^. 

We hrst present the estimates for the corrector x*" and briefly recall the proof. For detailed 
argument, we refer to m Lemma 4.1]. 


Lemma 4.1. Let d 
that 


2,3. Assume that [Al][A2][A3] and [S] hold. Then there exists C > 0 such 


nix^Wh < Ce'^- 


(4.3) 


Proof. The mean square of the L‘^{D) norm of the function is given by 

'E||X^lli2(D) = mGul^eU\\l2(^D) Id3 

Note that Ei'£{y)ne{z) = R^^iy — z) and use the bounds on the Green’s function. We get, for d = 3, 


EllX 


e\\2 


< 


L 


c 


Jd3 jic — y\'^~‘^\x — z\‘^~'^ 


R 


y- z 


\u{y)u{z)\dydz. 


In view of the bound (13^ . the product of Gu{-,y) and Gui-,z) is integrable on D with bounds 
independent of y or z. Hence, we can integrate over x first and then change variable in the remaining 
integral, which yields a factor of and verifies the desired result. The situation of d = 2 can be 
treated in the same manner. □ 


Next, we move on the remainder z^ = — u — x^- In view of the equations satisfied by u^, u 

and x'’) we have 

(-A + qe)z’^ + f{u^) - f{u) = -VeX^ + f{u)x^. 

Let /ig be the function 

f{u^) - f{u)^ 
u^-u 

where 1 denotes the indicator function. Then, in view of / € C^(M) and the bound (13.4p . we have 
he G L°°{D). We verify that z^ solves the Dirichlet problem 

- + qeZ^ + heZ^ = U^x" + {f{u) - he)x^, (4.4) 

with zero boundary condition. Indeed, this problem has unique solution and we verify easily that 
— u — x^ the solution. We have the following result. 

Lemma 4.2. Assume that the conditions in Lemma \4.1\ hold. Then there exists G > 0 such that 

mz^\\l2 < Ce^. (4.5) 
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Proof. Let denote the random linear differential operator — A + + hg). In view of [A3], the 

uniform bound on and the definition of h^, we observe that for each fixed e € (0,1) and w G 
the potential term := satisfies 


-Ai + c < / < M'. 

By standard elliptic theory, : Hq —)• is invertible. In particular the inverse operator 

(££,cj)-i jg g bouuded transformation on L?‘{D). In fact, < c~^. Applying 

this fact, we have 

Ik^|li 2 < 2c" 2 (||z^£X^lli 2 + Wifiu) - he)x^\h) ■ 


By [A2] and [A3], llz^gjjioo < C uniformly in e and iv. By the uniform bounds on u^,u, and by the 
smoothness of /, we confirm that \\f'{u) — h^\\L°° < C also uniformly in e and w. As a result, we 
have 11.2;^11^2 < C'llx^lli 2 - The desired result follows. □ 


The following control of tt'" — tt follows immediately. 

Corollary 4.3. Under the same conditions of Lemma \4.1\ there exists C > 0 such that 

-u\\l2 < Ce^. 


(4.6) 


5. Limiting Distributions of the Homogenization Error 

In this section, we identify the limiting distribution of the scaled homogenization error, i.e. the 
limit of the law of — u) in the space of Lp‘{D) functions. 

5.1. Expansion formnla and overview of proof. To characterize the limiting distribution of 
e “2 (yp — u), we extend to nonlinear equations the framework developed in PEI EE! for linear equa¬ 
tions. The first step is to find a suitable expansion formula similar to (|1.3p for the homogenization 
error. 

We observe that — u — satisfies the linear equation 

-Az" + gz" + f'iu)z^ = - (/(n") - /(u) - f'{u)e)- 

Hence z^ has the following representation 

= -QuT^ei^ - Gu {fiu^) - fiu) - f'{u)f,^) . (5.1) 

Substituting this formula to — u = x^ + and using the formula x^ = —Gut^eU, we obtain 
-u = - GuT^eU - Gul'eiu^ - u) - ^„(/(u^) - f{u) - f'{u)C^) 

— Gu^e^ T Gu^eGu^e'^ T Gu^eGu^eid^ ^) (^•^) 

- Guif{u^) - f{u) - f'{u)f,^) + GuVeGu{f{u^) - fiu) - /'(u)^^). 

With this expansion formula at hand, we will find the limiting distribution of £~^iu‘^ — u) by 
examining the terms on the right hand side one by one. In particular, we note that the first three 
terms are the same as those obtained for linear equations, while the last two terms involve the 
nonlinearity. 

We first recall the standard criterion for establishing weak convergence of probability measures 
determined by random processes that are L^(D) functions. 
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Theorem 5.1. Let {X^}, e G (0,1), he a family of random processes on the probability space 
(H, P) and C LP‘{D). Then converges in distribution in Lp‘{D) to the random process 

X in L^{D) if 

(i) (Convergence of inner products with test functions). For any (p G L?‘{D), the random 
variable {(p,X^) converges in distribution in M to {(p,X). 

(ii) (Tightness of distributions). The family of distributions in L'^{D) determined by {X^^} is 
tight. 

This Prohorov type result is well known and we refer to |23[ Chapter VI, Lemma 2.1] for a proof 
in the general Hilbert space setting. When the tightness of the distribution of {X^} is concerned, 
the following result becomes handy provided that {V^} in fact admits higher regularity. 

Lemma 5.2. Let {X^} be as in Theorem \5.1\ Suppose further that {X^} C H^{D), for some 
s G (0,1). Then the family of distributions in L‘^{D) determined by is tight if there exists 

some constant C > 0 such that 

E||X^|U»p) < C. (5.3) 

We refer to m Theorem A.l] for a proof. Since X^ is e 2 — u) in this paper, it indeed belongs 

to the more regular space H^{D). Nevertheless, to obtain the control (|5.3n . which is uniform in e, 
requires some work. 


5.2. Liming distribution for the homogenization error. The first three terms in the expan¬ 
sion formula (|5.2p are precisely those encountered in the linear setting; compare with (|1.3I) . We 
recall the following results. 


Lemma 5.3. Assume that the conditions of Theorem \2.1\ hold. Then as £ ^ 0, we have 

distribution ^ / \ / \ njr /r 

- 1 = - Gu{x,y)u{y)dWy, (5.4) 

ve“ JD 

with convergence in distribution in the space L‘^{D). Moreover, — u) converges in 

L^{FI,LP‘{D)) and Qu^eGu^eU converges in LP‘{FI,LP‘{D)) to the zero function. 


We briefly explain how these results are proved and refer to [181 Lemmas 4.2, 4.4, 4.5, 5.1 and 
5.3] for detailed proofs. 

(1) To show that —£~^Gui^eU has the correct limit, according to the criterion Theorem 15.11 
Remark 12.31 and Lemma 15.21 it suffices to show, for all ip G L?‘{D), 


^{-GuVeU,p) = 


Jd 


I X \ , , / \ , distribution 

n [—,uj ] u[x)m[x)dx - 

' e / £^o 


> Af (0, CJ^) , 


(5.5) 


where m := GuT and for some C > 0 and s G (0,1) independent of e. 


E 


e '^Gul^eU 


H‘ 


< C. 


(5.6) 


We recall that (j5.5p follows from a generalized central central limit theorem for oscillatory integrals 
with short range correlated random fields; see [H Theorem 3.7]. The uniform estimate (|5.6p for 
s can be found in [18] Lemma 5.1]; see Lemma 5.4 below where such a control is needed to estimate 
the first nonlinear term in (15.2|] . 
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(2) To show that Qu^eQu^e{u^ — 'u) converges to the zero function in {Q, L‘^(D)), we need the 
fact that for d = 2,3, 

nSul^eGuWh^L^ < Ce^ (5.7) 

where \\Gui'eGu\\L'^^L'^ is the operator norm. This can be proved by using the Green’s function 
bound p.7p and assumption [S]; see [TSl Lemma 4.5]. We then get 

WGul'eGul'eiy!' -'w)IIl2 < ||^«l^£^«||L2^i2||l/£||Lcx>||M^ -'u||l2. 


Take expectation, apply Holder inequality and then and Corollary 14.31 to get the desired 

estimate. 

(3) To show that Gut^eGut^eU converges to zero in L‘^{D)), the argument above is not valid 

because ||i 2 £m ||£2 is possibly of order one. To circumvent this lack of control, we note that 

WGul'sGul^eUWh < 2 {llGul^sGul^eU-^iGul^eGul^ey^Wh + \\E{Gul^eGul^eU)\\‘i2) ■ 

For the mean function ¥.{Gut^eGu’^e''J')i have 

W^iGu^eGuT^eu)]]^ = J ^Gu{x,y)Gu{y, z)R u{z)dzdy^ dx 


■- [ Gu{x,y)Gu{x,i)Gu{y,z)Gu{i,y)R 
Jd^ 


y- 


u{z)u{r])d^dr]dydzdx. 


Integrate over x first and note that, for d = 3, ||u||loo < G and Gu is square integrable and 

1 


'D 


\Gu{x,y)Gu{x,0\dx<G [ 

Jl 


d\x- y\^ “^Ix - 


dx < G. 


d-2 — 


We get 


\\^{Gut^eGul^eU)\\l2 <g[ 

JD- 


1 


R 


y - z\ 


d^dydydz. 


Id^ \y - ^ 

After a change of variables and using the fact that 7?(-)/| ' is integrable over W^, we verify that 
the above integral is of order <C e^. It follows that \\s~2K{GuI^£GuI^£u)\\l^ converges to zero. The 
case of d = 2 is similar. 

For the variation I 2 := Gu^eGu^eU — ^.{Gu^eGu^eu), we have for d = 3, 

2 


E||/ 2 |li 2 =E 


D \JD^ 


Gu{x, y)Gu{y, z) [Ve{y)Ve{z) - "&Ve{y)ve{z)] u{z)dzdy dx 


[ Gu{x, y)Gu{x, y')Gu{y, z)Gu{y', z')u{z)u{z')^ 
Jd^ 


,^l,^XX]dz'dy'dzdydx. 
see e 


in 


Now, using the estimates for given by (|2.8p . and by standard techniques for potential inte¬ 
grals, we get E||/ 2||£2 < Ge^'^ <C e'’*, which shows that e~^l 2 converges to the zero function i 
L2(H,L2(T))). For d = 2, we have the same result. 

Now we deal with the two nonlinear terms in (|5.2p . 


Lemma 5.4. Under the assumptions of Theorem \2.1{ as e ^ 0, we have 

Gu{f{u^) ~ f{'^) ~ — u)) distribution^ ^ 

V? 

with eonvergence in distribution in L?‘{D) and 0 denotes the constant zero funetion. Moreover, the 
term e~^Gu^^eGuif {u‘^) — f{u) — f'{u){u‘^ — u)) converges to the zero function in L^{yt.,L?‘{D)). 
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Proof. Part I: The convergence of QuVeGu{f{u^) — fiu) — f'iu)iu^ — u))- This part is simpler and 
is similar to the control of the remainder term QuT^eGu^e{u^ ~ u) in the linear setting. 

Firstly, by Taylor’s theorem, the regularity of / and the uniform (in e) bound of and u, we 
have 

l/(^'') - f{u) - f'{u){u^ - u)| < C'\u^ - n|^, (5.9) 

where C = ||/||c 2 ([- 2 C, 2 C]) C is the bound in p.4l) . This shows 

- f{u) - f'iu){u^ - u)\\l2 < C\\u^ - u\\l2. 

Applying the uniform bounds on the operator norm of Gu^'eGu, we have again 

WGut^eGuif (u^) - fiu) - /'(n)(n^ - u))\\l 2 < CWGuJ^sGuWl^^l^Wu^ - u\\l 2 . 


The desired result then follows by taking expectation and applying the Holder inequality, the bound 
(15.71) and (14.6p . 

Part II: The limiting distribution of e~ 2 Q^(^f{u'^) — f[u) — f'{u){u^ — u)). We follow again the 
criterion in Theorem O 

Convergence in distribution of inner products. Fix an arbitrary ip G L^iD). We need to consider 
the distribution of 2 (^„(/(ti^) — f(^u) — f'{u)ff,ip) = e~ 2 {f{u^) — f{u) — f'{u)ff ,m), where 
'm = QuT- Note that ||m||x,oo < C'||y5||^2. Using (15.91) . we get 

|(m,/(n^) - fiu) - f'iu)iu^ - u))| < C\\ip\\L 2 \\u^ - u\\\ 2 . 

As a result, we get 


E 


e 2 (m,/(u^) —/(u) —/'(u)(u'^ — u)) < Ce 2 E||u^ — u ||^2 < C'e 2 . 


(5.10) 


This show that (e^Quifiu^) —fiu) —f'iu)iu^— u)),(p) converges, inL^(H) and hence in distribution, 
to 0, agreeing with the trivial distribution of (v^,0). 

Tightness. Let = /(u^) — fiu) — f'iu)f^^. The function we are considering is QuC. To show 
tightness of the distributions of {Guf^}, we control the expectation of the norm of Gu^^ for some 
s G (0,1). For the semi-norm, we calculate 

2 1 f \Gur%x) -GuCiy)\‘^ 


£ ^GuT^ 


J H‘ Jd2 


_7/1 d-\-2s 


\x - y\ 


-dydx 


1 f iGuix,z) - Guiy,z)) iGuix,r]) - Guiy,v)) 


£‘^ JD'i 


r^iz)r^ir]) dzdydydx. 


\x — 

It is proved in m Lemma 5.1] that the uniform bounds (|3.7p on the Green’s function and its 
gradient imply that 

\iGuix,z) - Guiy,z))iGuix,r]) - G„(y,r/))| 

Id^ 


L 


-dydx < C, 


(5.11) 


\x — 

for any s G (0,^), uniformly for z,r] G D. We note also that (15.9p implies, for some C > 0 
independent of e and w, 

llr^llLi < C\\u^ - n|||2. 

Now we integrate over x and y first and use the inequalities above, and get 


E 


e "^GuT^ 


12 


J H= 


< Ge-^^ ( / yiz)\dz ) < Ge-^ ||r^||LooE||r^||ii < G. 


'D 


—d ||„£ I 
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Similarly, one can show that K\\s~^Gur‘^\\L^ < C as well. We conclude that {K\\s~^Qur^\\H‘’{D)} is 
uniformly bounded. It follows that the family of probability distributions in L^{D) determined by 
{GuT'^} is tight. We get the desired result by an application of Theorem 15. II □ 

Finally, we combine all of the results above and prove the second main theorem. 

Proof of Theorem \2.2[ According to the above results, in the expansion formula (j5.2|] . the first term 
converges in distribution in the space of L?‘{D) functions to the desired limit of Theorem 12.21 All 
other terms converge in distribution to the zero function, which is deterministic. Therefore, these 
terms converge to zero also in probability. As a result, the sum of all terms converge in distribution 
to the limit of the leading term. This concludes our proof. □ 

6. Further Discussions 

6.1. Long range correlated random fields. We have assumed so far that the random potential 
q{x,uj) had short-range correlation. Our approach applies to the setting of some long range cor¬ 
related potentials as well. Following [HE], a large class of long-range correlated potential can be 
constructed as functionals of long range correlated Gaussian random fields. Let q{x, oj) = q+i'{x, oj) 
with q a nonnegative constant; we assume 

[LI] u(x, uj) = $( 0 (x)); g(x, uj) is a centered stationary Gaussian random field with unit variance. 
Furthermore, the correlation function of 0 (x,a;) has a heavy tail in the sense that, for some 
positive constant Kg and some real number a E (0, d), 

Rq{x) := E{ 0 (i/, io)Q{y -|- x, w)} ~ Kg|x|“" as |x| —>■ oo. ( 6 . 1 ) 

[L2] The function $ : M ^ M satisfies —M < ^ + q < M, and has Hermite rank one, i.e. 

/ <h(s)e“~fis = 0 , Vi := / 7 ^ 0 . ( 6 . 2 ) 

Jr. Jr 

[L3] The Fourier transform of the function <1> satisfies 

[ 16(01(1 + 1^1') < 00 . (6.3) 

Jm 

We henceforth refer assumption [L] to the above conditions all together. 

The assumption [L2] ensures that u(x,a;) = $( 0 (x,a;)) is mean zero and the bounds on <1> ensure 
that |g(x,a;)| < M, which is (12.2p . From the above construction, we check that u(x,w) is stationary 
ergodic and has a long-range correlation function decaying as k|x|“", k = V^Kq] see [THl Lemma 
A.3] for the details. Assumption [L3] allows one to derive a (non-asymptotic) estimate, see [THl 
Lemma A.5], for the fourth-order moments of v{x,oj). We have the following analog of Theorem 
12.11 and Theorem 12.21 

Theorem 6.1. Let d = 2,3, and u be the solutions to dEU) and mi respectively. Suppose that 
[Al][A2][A3] and [L] hold, g E L^[D). Then, there exists positive universal constant C, such that 

E\\u^ - u\\l 2 < C{\\g\\L2)e^. (6.4) 

Let W^{dy) be defined formally as W°‘{y)dy and W°‘{y) be a centered stationary Gaussian 
random field with covariance function k\x — y\~°‘, where k = KqVi >0 and Kg and VT are defined 
as in dS+D and (j6.2p . Then we have the following. 
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Theorem 6.2. Suppose that the assumptions in Theorem \ 6.1\ hold. Let Gu{x,y) he the Green’s 
function of (j3.6p . Then as e ^ 0, 

^ distribution^ ^ r y)u{y)W-{dy), ^n L\D). (6.5) 

Ve" Jd 

Remark 6.3. We will not present the proof since they can easily be adapted from the approach 
in mm and our control of the nonlinear terms earlier. The right hand side of (16.5p is an integral 
with respect to the multiparameter Gaussian random processes W^, we refer to |19] for the theory. 
Let us denote this integral by X] then it determines a Gaussian distribution in the space L?‘{D). 
In particular, for any {p G the inner product {ip,X) has Gaussian distribution 

with 

2 f iuGuT)iy)iuGuT)iz),, 


6.2. Non-separated nonlinearities. We have so far assumed that the heterogeneous reaction 
function was of the form f^{x,u,uj) = q£{x,uj)u + f{u)] in other words, the random potential 
and the nonlinear function / are separated. This choice is made mainly for notational simplicity. 
By some careful modifications of the main assumptions, we can extend our result to a class of 
general nonlinearities in non-separated form and consider a heterogeneous problem of the form: 

( - Au^ + f = g{x), xGD, 

I = 0, X G dD. 

Here, the nonlinear reaction is /^(x, u, to) = f{j,u, u). We may next modify [Al] [A2] [A3] as follows. 
[Al’j f{x,s,Lj), X G s G K, is a random field on P) and stationary ergodic in x. In 

particular, there exists f{s,co) such that f{x,s,uj) = f{s,Tx 0 j). 

[A2’] There exist M > 1 and 7 > 1 such that 

|/(s,t<;)| < M(1 -|- jsp), for all s G M and cu G H. ( 6 . 8 ) 

Moreover, if d = 3, we further assume that 7 < 2d/{d — 2) — 1 = 5. 

[A3’] f{-,uj) : M —)• M is continuously twice differentiable uniformly in u, and for some c > 0 , 
71,72 > 0, for all a; G H and s G M, 

Al + 7'(s,w) > c, (6.9) 

l7'(s,a;)l < M(1 + Js]^^), \fis,u;)\ < M(1 + Js^^). (6.10) 

Here f'{s,u) and f''{s,u}) denote, respectively, the first and second order derivatives of / with 
respect to s. Let f{s) be the mean E/(y,s,-) = E/(s, •). Then the mean-zero fluctuation of the 
random potential, i.e. n{y,Lo) defined in (12.4p . is replaced by 

’^iy,s,uj) = f’{y,s,uj)-J'{s). ( 6 . 11 ) 

We also let fi{y,s,u}) := f{y,s,Lj) — f{s), which is the antiderivative of n and fj,£{x,u) plays the 
role of i's{x)u in the separated case. 

To ensure that the random medium has short range correlations, we assume that: 

[S’] There exists g : M+ —M satisfying [S] such that, for each s G M, the maximal correlation 
function of the family of random fields {i^(x, s), /i(x, s)}, x G K, is bounded above by g. 
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Under [Al’][A2’][A3’] and [S’], the approach presented in this paper may be applied to prove 
that of (16.7p converges weakly in and strongly in as e ^ 0 for a.e. a; G fl, to the solution 
u of 


-/^u +f{u) = g{x), xeD, 

u = 0, X G dD. 


( 6 . 12 ) 


Indeed, the existence and well-posedness of the above equations are guaranteed by [A2’][A3’]. We 
note that the proofs of Lemma 13.11 and Lemma 13.21 can be adapted since they only require the 
uniform lower bound on the derivative of / and the uniform bound on growth rate of /. Hence, 
we still have uniform L°° bound on . The non-separated form does introduce some technicality 
in the proof of almost sure homogenization theory because, in the application of ergodic theorem, 
we have to deal with the u-dependent random process f{y,u,u}). Even though converges 

weakly in for all p G [l,oo) to /(s) a.e. in U for each fixed s, to show /(|,m(x)) converges 
to f{u{x)), one needs extra effort. For instance, a full measure event fli can be constructed by 
intersecting over rationals s G Q, where is the event on which /(f , s) converges to /(s). The 
proof of Lemma 13.31 then gets through by an additional approximation of u using simple functions 
with rational values. 

Once homogenization theory is established, we can continue to carry out the error estimates and 
the analysis of the distribution of fluctuations. The linearized equation, i.e. the analog to (13.6p . is 


- AGu{x; y) + /iu{x)) Gu{x] y) = 6y, for x G D, 

Gu{x]y) = 0, for x G dD. 


(6.13) 


The leading term in the distribution will be given by the solution to 

+l\u)x^ = -Pe{x,u,uj), in H, 


which can be conveniently written as = —QulJ-eix, u). Here, Qu denotes the fundamental solution 
operator to (j6.13p . We then get the expansion formula 

-U = -Qu {p,eix,u)) - QuV^{x,u){u^ - u) - Qu (/e(x,U^) - fe{x,u) - f'^{x,u){u^ 

= -Qu {iJ-e{x,u)) + Qul'eix,u)QuPeix,u) + GuJ^e{x,u)Qul'e{x,u){u^ - u) 

- Qu {feix, U^) - feix, u) - f'^{x, u)^^) 

+ QuT^e{x, u)Qu {fe{x, U^) “ fe{x, u) - f^{x, u)^^) . 

In view of the modified assumption [S’], and by following the argument developed in earlier sections, 
we show that all items on the right except the first one are of negligible order in distribution. In 
particular, as an analog to Theorem 12.21 we have 

distribution^ / Gu{x,y)a,{u{y))dW{y), in L\D). (6.15) 

Ve“ JD 

Here, cr^(s), for each fixed s G M, is defined as 

(T^(s) = / Ry{x,s)dx= / IE/i(x, s)/i(0, s)(ix. 

JR'* 


-u)) 

(6.14) 
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6.3. Further studies. We conclude this section by mentioning a couple of extensions to the above 
results. First, this paper considers only the physical dimension d = 2,3, and it would seem natural 
to extend the studies of this paper to arbitrary dimension d > 4. For linear equations, as studied 
by [18] , the framework of mum can still be applied, and is more or less unchanged, provided that 
we seek for the limiting distribution in or other weaker spaces. This approach may not apply 
directly in the nonlinear setting. The uniform estimates on is essential in our treatment of 
the nonlinearity, and such an estimate is not available in higher dimension for general source term 
g, the right hand side in 11.11 in and not more regular. It remains to explore how to generalize 
the framework to other functional settings and to develop new ways to control the nonlinearity 
terms. 

Another related further study is to modify the analysis in the continuum setting to handle 
the discrete setting, i.e., to address the numerical methods of m- As shown in [SIEI for the 
linear equations, studying the limiting distribution of the difference between the solutions to the 
heterogeneous equation and to the homogenized solution, obtained from multi-scale numerical 
schemes, in the limit of e ^ 0 and then the discretization size h ^ 0, and comparing the results 
with the theory in the continuum setting, one can build a benchmark to assess the performance 
of numerical schemes in capturing numerically the fluctuations of heterogeneous equations. To 
perform such an analysis in the nonlinear setting requires new ideas, in addition to those in [SIE]. 


Acknowledgements 

GB acknowledges partial support from NSF grant DMS-1408867. WJ acknowledges partial 
support from NSF grant DMS-1515150. The authors thank the referees for their helpful comments. 


References 

[1] G. Bal. Central limits and homogenization in random media. Multiscale Model. SimuL, 7(2):677-702, 2008. 

[2] G. Bal, J. Gamier, Y. Gu, and W. Jing. Corrector theory for elliptic equations with long-range correlated random 
potential. Asymptot. Anal., 77(3-4):123-145, 2012. 

[3] G. Bal, J. Gamier, S. Motsch, and V. Perrier. Random integrals and correctors in homogenization. Asymptot. 
Anal, 59(l-2):l-26, 2008. 

[4] G. Bal and W. Jing. Homogenization and corrector theory for linear transport in random media. Discrete Contin. 
Dyn. Syst, 28(4): 1311-1343, 2010. 

[5] G. Bal and W. Jing. Corrector theory for elliptic equations in random media with singular Green’s function. 
Application to random boundaries. Commun. Math. Sci., 19(2):383-411, 2011. 

[6] G. Bal and W. Jing. Corrector theory for MSFEM and HMM in random media. Multiscale Model. SimuL, 
9(4): 1549-1587, 2011. 

[7] G. Bal and W. Jing. Corrector analysis of a heterogeneous multi-scale scheme for elliptic equations with random 
potential. ESAIM Math. Model. Numer. Anal., 48(2):387-409, 2014. 

[8] G. Bal and K. Ren. Transport-based imaging in random media. SIAM J. Appl. Math., 68(6):1738-1762, 2008. 

[9] A. Bourgeat and A. Piatnitski. Estimates in probability of the residual between the random and the homogenized 
solutions of one-dimensional second-order operator. Asymptot. Anal, 21(3-4):303-315, 1999. 

[10] T. Cazenave. An introduction to semilinear elliptic equations. Editora do IM-UFRJ, Rio de Janeiro, 2006. 

[11] E. Di Nezza, G. Palatucci, and E. Valdinoci. Hitchhiker’s guide to the fractional Sobolev spaces. Bull. Sci. Math., 
136(5):521-573, 2012. 

[12] R. Figari, E. Orlandi, and G. Papanicolaou. Mean field and Gaussian approximation for partial differential 
equations with random coefficients. SIAM J. Appl. Math., 42(5):1069-1077, 1982. 

18 


[13] A. Gloria, S. Neukamm, and F. Otto. Quantification of ergodicity in stochastic homogenization: optimal bounds 
via spectral gap on Glauber dynamics. Invent. Math., 199(2):455-515, 2015. 

[14] A. Gloria and F. Otto. Quantitative results on the corrector equation in stochastic homogenization. Preprint, 
2014. 

[15] Y. Gu and J.-C. Mourrat. Scaling limit of flucutations in stochastic homogenization. Preprint, 2015. 

[16] M. Hairer, E. Pardoux, and A. Piatnitski. Random homogenisation of a highly oscillatory singular potential. 
Stock. Partial Differ. Equ. Anal. Comput., l(4):571-605, 2013. 

[17] V. V. Jikov, S. M. Kozlov, and O. A. Oleinik. Homogenization of differential operators and integral funetionals. 
Springer-Verlag, Berlin, 1994. 

[18] W. Jing. Limiting distribution of elliptic homogenization error with periodic diffusion and random potential 
(long version). arXiv:1505.02721, 2015. 

[19] D. Khoshnevisan. Multiparameter proeesses. Springer Monographs in Mathematics. Springer-Verlag, New York, 
2002. An introduction to random fields. 

[20] J.-C. Mourrat and J. Nolen. Scaling limit of the corrector in stochastic homogenization. Preprint, 2015. 

[21] J.-C. Mourrat and F. Otto. Correlation structure of the corrector in stochastic homogenization. Preprint, 2014. 

[22] J. Nolen and G. Papanicolaou. Fine scale uncertainty in parameter estimation for elliptic equations. Inverse 
Problems, 25(11):115021, 22, 2009. 

[23] K. R. Parthasarathy. Probability measures on metric spaces. Probability and Mathematical Statistics, No. 3. 
Academic Press, Inc., New York-London, 1967. 

Department of Applied Physics and Applied Mathematics, Columbia University, 500 W. 120th Street, 
NY 10025, USA 

E-mail address: gb2030@coluinbia.edu 

Department of Mathematics, The University of Chicago, 5734 S. University Avenue Chicago, IL 
60637, USA 

E-mail address: wjing@math.uchicago.edu 


19 



